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PREFACE 

This  report  is  a  replacement  for  the  report,  "The 
Dynamic  Design-Analysis  Method  (DDAM)  in  NASTRAN,"  David  W. 
Taylor  Naval  Ship  Research  and  Development  Center  Report 
DTNSRDC-81/073,  October  1981.  The  present  work  incorporates 
efficiency  improvements,  corrections,  and  additional  capa¬ 
bility  into  the  previous  work,  necessitating  modifications 
to  the  user  specifications.  This  report  is  self-contained 
and  does  not  require  retention  of  DTNSRDC-81/073. 
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ABSTRACT 


This  report  describes  the  theory,  implementation,  and  use 
of  the  Dynamic  Design-Analysis  Method  (DDAM)  in  the  NASTRAN 
finite  element  structural  analysis  computer  program.  DDAM  is 
the  procedure  used  in  the  shock  design  of  shipboard  equipment. 

Since  such  equipment  is  also  frequently  analysed  with  NASTRAN, 
the  inclusion  of  DDAM  in  NASTRAN  greatly  enhances  the  effi¬ 
ciency  of  the  design-analysis  process. 

ADMINISTRATIVE  INFORMATION 

The  work  reported  was  performed  under  Tauk  Area  S0923001,  Task  23350,  Work  Unit 
1102-020. 

INTRODUCTION 

The  Dynamic  Design-Analysis  Method  (DDAM)  is  the  standard  procedure  for  shock 

design  of  shipboard  equipment.  Often,  the  equipment  is  first  analyzed  with  the 
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NASTRAN  finite  element  structural  analysis  computer  program.  The  data  and  results 
must  then  be  converted  into  other  forms  for  use  in  DDAM.  Incorporating  DDAM  into 
NASTRAN  has  enabled  the  entire  process  to  be  performed  more  efficiently.  This 
report  describes  Che  implementation  of  DDAM  into  the  DTNSRBC  version  of  NASTRAN  Level 
17.5  and  the  input  data  required . 

The  present  work  was  undertaken  as  part  of  an  effort  to  develop  validated 
structural  analysis  procedures  in  support  of  the  Improved  Performance  Machinery 
Program,  a  DTNSRDC  Systems  Development  Department  project  related  to  onboard 
machinery, 

THEORY 

The  steps  of  the  DDAM  procedure  are  described  here  very  briefly.  A  more 
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complete  description  is  gLven  by  Belshelo  and  O'Hara.**  Their  terminology  as  well  as 
that  of  NASTRAN  will  be  used  in  the  mathematical  description. 

Step  1.  Compute  natural  frequencies  and  mode  shapes.  (NASTRAN's  Rigid 
Format  3) . 

Step  2.  Compute  the  participation  factor  for  each  mode. 


*A  complete  listing  of  references  is  given  on  page  19. 
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where  is  the  mass  of  the  ith  degree-of-freedom  and  X^a  Is  the  ith  component  of 
the  ath  mode  shape. 

It  is  assumed  that  only  those  terms  of  {X  }  that  correspond  to  a  particular 
direction  are  used  in  Equation  (1).  That  is,  the  ath  mode  may  have  three  participa¬ 
tion  factors  associated  with  it,  one  for  each  orthogonal  direction. 

In  NASTRAN  terminology,  the  numerator  of  Equation  (1)  may  be  written  as 
(considering  all  computed  modes) 


[*]T  [M]  [V]  (la) 

where  [$]  *  matrix  of  eigenvectors  (mode  shapes),  order  n  x  m,  with  n  =■  order  of  the 
problem,  m  ■  number  of  modes  computed 

(M]  *  mass  matrix,  order  n  x  n 

[V]  ■  direction  cosine  matrix,  order  n  *  l  with  l  **  1,  2,  or  3,  the  number 
of  desired  directions 

Typically,  (V)  may  consist  of  lrs  and  0's  which  "piek  off"  desired  directions. 
However,  that  form  is  not  nocessary  and  any  consistent  set  of  direction  cosines  may 
be  used. 

The  denominator  of  Equation  (1)  may  be  written  as  (considering  all  computed 
modes) 


($]T  in)  m  (ib) 

which  is  the  diagonal  modal  mass  matrix.  Therefore,  we  may  write  Equation  (1)  as 
(PaW  *  Idla*  <WXWH*)>rl  l*)T  (N)  (V)  (2) 
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Step  3.  Compute  the  "effective  mass"  in  each  mode 


M 

a 


P  Z  M.  X. 
a  i  i  ia 


(3) 


In  NASTRAN  terminology, 

>\W  =  >V  *  [W1  I"1  EVl]  «> 

where  the  x  on  the  right  side  indicates  the  so-called  matrix  outer  product,  in  which 
a  term-by-term  multiplication  is  performed.  For  example,  if 

[Cl  -  (A)  x  [B] 


then 
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Step  4.  Compute  the  "effective  weight"  [ W_^ ]  in  each  mode  by  multiplying  ] 
by  g,  the  acceleration  due  to  gravity. 

Step  5.  Compute  the  direction-dependent  velocity  spectrum  design  values  (V  ) 

from  (W  ]. 
a 

Step  6.  Compute  the  effective  static  load  at  each  mass,  due  to  the  ath  mode ,  bv 


M, 


V 


a 


V  co 

U  *1 


(5) 


where  cu  is  the  ath  natural  radian  frequency. 

u 

In  NASTRAN,  the  matrix  of  loads  is  computed  as  follows:  The  matrix  product 
(M)l$)  Is  of  order  n  *  m  and  corresponds  to  the  product  M^X^  of  Equation  (5). 
(Here,  only  terms  of  in  the  desired  directions  are  used.)  The  ath  column  of 

corresponds  to  the  ath  mode.  Multiplying  the  ath  column  by  «  and  by  P  V 

ii  M 
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for  the  first  desired  direction  gives  a  matrix  of  load  vectors  of  order  n  x  m.  If 

P  V  *s  for  other  desired  directions  are  used,  then  other  n  *  m  sets  of  loads  are 
a  a 

created  and  appended  to  the  first  set.  A  final  load  matrix  [F],  of  order  n  x  mfl,,  is 
thus  created,  where  i  is  the  number  of  desired  directions,  i.e.,  there  will  be  mi 
static  load  cases.  (In  practice,  instead  of  the  product  V  to  in  Equation  (5),  the 

3  3 

term  actually  used  is  min  (V  to  ,A  g),  where  A  is  an  acceleration  spectrum  design 
value  in  g's  and  g  is  the  acceleration  due  to  gravity.) 

Step  7.  Perform  static  analyses  to  compute  direction-dependent  maximum 

responses,  using  the  load  cases  from  Step  6,  and  calculate  element 
stresses. 

The  computation  of  the  effective  static  load  at  each  mass  in  Step  6  and  the 
static  analyses  of  Step  7  may  be  replaced  as  follows:  For  the  ath  mode  {$  }, 
Equation  (5)  may  be  written 


[M]  (*„)  Pa 


V 

a 


w 

a 


•  (5a) 


Solving 


(K]  {u  )  -  {F  }  (5b) 

a  a 

where  (K|  Is  tho  stiffness  matrix  and  {u  }  is  the  vector  of  grid  point  displacements, 

a 

yields 

(u  )  -  IK)-1  (H)  {*  )  «  (Sc) 

However,  from  dynamics, 

t-*J l  M+K]  {♦  )  -  0  (5d) 

a  3 

or 

U)“l  CHS  ($a)  -  ^  ($a)  (5e) 
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Using  Equation  (5e)  in  Equation  (5c)  yields 


P  V 


a 


(5f) 


(As  in  Step  6,  rather  than  V  / uj  ,  min  (V  /u>  ,A  g/w  )  is  used.) 

a  a  a  a  a  a 

Equation  (5f)  is  used  in  NASTRAN  to  compute  the  direction-dependent  maximum 
response. 

Step  8.  For  each  of  the  £  desired  directions,  compute  the  NRL  sums  of 
3 

stresses  for  each  element  as  follows: 


S 


j 


where  ■  maximum  stress  at  the  Jth  point  (taken  over  the  modes  under  consideration 

in  one  desired  direction)  and  ■  stresses  (other  than  the  maximum)  at  the  jth 

point  corresponding  to  the  modes  described  for  S.  . 

jo 

NASTRAN  IMPLEMENTATION 

Since  DDAM  requires  the  determination  of  natural  frequencies  and  mode  shapes, 
a  NASTRAN/ DDAM  analysis  will  consist  of  Rigid  Format  3  with  ALTER’ s.  The  ALTER*  si 
required  in  the  NASTRAN  run  will  compute  the  various  quantities  described  in  the 
previous  section.  Among  the  ALTER’s  are  instructions  for  NASTRAN  to  execute  the 
following  eight  new  functional  modules  which  were  added  to  the  program  for  DDAM. 


CENCOS 

GEKCOS  generates  the  direction  cosine  matrix  (V)  in  Expression  (la).  The  user 
may  specify  a  coordinate  system  which  defines  tie  shock  directions.  A  PA RAM  bulk 
data  card  giving  the  value  of  parameter  SHOCK  passes  to  GEKCOS  the  coordinate  system 
identification  number  of  the  shock  system.  If  the  user  does  not  include  such  a  card, 
HASTRAN's  basic  coordinate  system  will  be  used.  (The  value  of  parameter  SHOCK 
should.  In  most  cases,  correspond  to  the  displacement  coordinate  system 
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identification  number  for  the  grid  points  in  the  problem.  However,  to  allow  for 
possible  exceptions,  no  check  for  this  correspondence  is  made.)  Parameter  DIRECT 
must  also  be  specified,  defining  the  directions  of  the  shock  system  which  are  to  be 
considered.  The  options  for  DIRECT  are  1,  2,  3,  12,  13,  23,  and  123.  ror  example, 
if  DIRECT  is  23,  then  the  second  and  third  directions  of  the  shock  coordinate  system 
will  be  used.  If  the  user  does  not  define  DIRECT,  the  default  is  123,  i.e.,  nil 
three  directions  will  be  considered. 

The  DMAP  statement  for  GENCOS  is 

GENCOS  BGPDT,  CSTM/DIRCOS/C.Y,  SHOCK»0/C,Y,  DIRECT-123/ 

V,N,  LUSET/V.N.NSCALE  $ 


DDAMAT 

DDAHAT  calculates  a  matrix  outer  product  such  as  that  in  Equation  (4),  and 
multiplies  the  result  by  parameter  GG.  For  example,  to  compute  effective  weights. 
Steps  3  and  4  are  performed,  and  GG  •  g  *  386.4  if  units  of  pounds  and  inches  arc 
used. 

The  DMAP  statement  for  DDAMAT  is 

DDAHAT  A,B/C/C,Y,CG  $ 

Parameter  GG  must  be  given  a  value  on  a  PARAM  bulk  data  care  or  in  the  DMAP  statement 
itself. 


CENPAKT 

It  is  assumed  that,  in  the  eigenvalue  analysis,  the  lowest  N  nodes  were  computed. 
If,  in  the  static  analyses  (or  equivalent  static  analyses),  fewer  modes  are  to  be 
used,  say,  the  lowest  M,  where  N  <  N,  then  the  orders  of  a  number  of  matrices  must 
be  truncated.  CENPART  generates  the  partitioning  vectors  which  KASTRAK  functional 
module  PARTH  requires  to  partition  the  necessary  matrices. 

The  DMAP  statement  for  GEKPART  is 

GENPART  PP/'RPLAHB,  CP  LAMB,  RPPF,CPMP/C,Y,mODES/V,JJ, RHODES  $ 
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Parameter  LMODES  is  the  integer  value  of  the  number  of  lowest  modes  to  be  used  in 
the  static  analyses.  The  value  of  this  parameter  must  be  specified  on  a  PARAM  bulk 
data  card,  or  else  a  NASTRAN  fatal  error  will  result. 

DESVEL 

DESVEL  computes  design  velocity  and  acceleration  spectra.  The  assumed  form  for 
velocity  is 


V  = 


V.  +  w 

V  - 

a  V  +  W 
c 


(6) 


where  V  =  velocity  computed  from  modified  effective  weight  W 
=  factor  usually  associated  with  a  desired  direction 
Va>V^,Vc  =  factors  usually  associated  with  various  ship  types  and  parameters 
W  =  effective  weight/1000. 

Items  V  and  are  in  units  of  length/ second,  and  and  Vc  are  in  units  of  effec¬ 
tive  weight  W. 

Acceleration  spectrum  values  may  be  expressed  in  one  of  two  forms.  The  first 
form  is  the  same  as  that  for  velocity 


A 


A.  +  W 

A  A  _ 

f  a  A  +  W 
c 


(7) 


The  second  form  is 


A  =  AfAa  (Ab+W)  (Ac+W)/(Ad+W)2 


(8) 


wh^re  A  =  acceleration  in  g's  for  modified  effective  weight  W.  Af,  A^,  A^,  Ac,  and 
A^  ere  factors  defined  similarly  to  factors  V^,  Vfl,  V^,  and  V£.  If  A^  =  0,  then  the 
form  of  Equation  (7)  is  used.  In  addition,  values  Vto/g  are  computed  and  are  output 
along  with  acceleration  values  A  for  comparison  purposes.  Also,  a  matrix  of  minimum 


values  of  Vw  and  Ag  is  output  for  use  in  Equation  (5),  i.e., 


Amin  =  min  (Vw‘AS) 


(9) 


Finally,  the  matrix 


A'  =  —  A  . 
min  ,  i  min 
<u 


(9a) 


is  output  for  use  in  Equation  (b£).  Note  that  the  natural  frequency  must  not  be 
zero.  However,  this  is  not  a  restriction  for  DDAM  since  a  fixed  base  is  assumed. 

The  DMAP  statement  for  DESVEL  is 

DES1.  t  E^FW , OMEGA/ SSDV , ACC , VWG ,MINAC,MIN0W2/C, Y , GG  =  386.4/ 
C,Y,VEL1/C,Y,VKL2/C,Y,VEL3/C,Y,VELA/C,Y,VELB/ 
C,Y,VELC/C,Y,ACC1/C,Y, ACC2/C,Y,ACC3/C,Y, ACCA/ 

C,Y, ACCB/C,Y,ACCC/C,Y,ACCD  $ 

Parameter  GG  is  the  acceleration  due  to  gravity.  A  default  value  of  386.4  is 
supplied.  Any  other  value  must  be  specified  on  a  PARAM  bulk  data  card.  Parameters 
VEL1,  VEL2,  and  VEL3  correspond  to  factor  in  Equation  (6)  in  the  first,  second, 

and  third  desired  directions,  respectively.  If  fewer  than  three  directions  are 
desired,  then  only  VELi,  or  VEL1  and  VEL2,  are  specified.  For  example,  if  only  one 
direction  is  specified,  say  direction  3,  then  VELI  corresponds  to  direction  3,  the 
first  (and  only)  desired  direction.  Parameters  VELA,  VELB,  and  VELC  correspond  to 
V  ,  V^,  and  V  respectively,  in  Equation  (6).  These  velocity  parameters,  VELI 
through  VELC,  must  appear  on  PASAri  bulk  data  cards,  or  else  a  NASTRAN  fatal  error 
will  result.  If  VEL2  or  VEL3  is  not  used,  then  values  of  0.  must  be  specified. 

Acceleration  parameters  ACC1  through  ACCD  are  similar  to  VELI  through  VELC  and 
refer  to  Equations  (7)  and  (8). 
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DDAMPG 

DDAMPG  creates  the  static  load  vectors  of  Equation  (5)  or  the  maximum  responses 
of  Equation  (5f ) .  For  the  former,  the  matrix  [KP]  *  [MJ[$]  is  input  to  DDAMPG  and 
is  operated  on  by  a  matrix 


[PVW]  -  [Pal  x  [Acinl 


where  [P^]  is  the  matrix  of  participation  factors  defined  by  Equation  (2)  and  [A  ^Q] 
is  computed  from  Equation  (9) .  The  order  of  these  matrices  is  m  *  l  where  m  io  the 
number  of  modes  to  be  used  and  i  is  the  number  of  desired  directions.  [PVW]  is 
formed  by  functional  module  DDAMAT. 

The  columns  of  [PVW]  correspond  to  desired  directions.  Within  a  column,  each 
row  term  corresponds  to  a  mode.  The  matrix  [HP]  is  of  order  n  *  m,  where  n  is  the 
number  of  degrees-of-freedom  in  the  problem.  Each  column  of  IMP]  corresponds  to  a 
mode,  and  in  each  column  of  [PVW],  the  ith  row  term  of  [PVW]  multiplies  the  ith 
column  of  [HP].  After  all  columns  of  [PVW]  have  been  considered,  the  resulting 
static  load  matrix  is  of  order  n  *  (mi). 

To  compute  the  maximum  response  of  Equation  (5f),  the  same  operations  just 
described  are  performed,  except  that  matrix  [PHIG]  *  [$]  replaces  [MP]  and 


IPVOW]  -  -j  lAnlnl 

(ii 


(10a) 


replaces  [PVW],  where  «  “  to  for  the  ath  row  of  [A  .  ]. 

a  mm 

The  DMAP  statement  for  DDAMPG  for  static  loads  is 


DDAMPG  MP,PVW/PG/V,N,NXODES/V,N,NDI&  § 


For  maximum  responses,  the  DMAP  statement  is 


DDAMPG  PHIG , PVOW,' UGV/V , N , NMODES/V , H , NDIR  $ 
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W3S TO 


CASEGEN 


The  static  load  and  maximum  response  vectors  created  in  DDAMPG  are  considered 
individual  load  cases  by  NASTRAN  and  must,  therefore,  be  selected  in  the  Case  Control 
Deck.  The  number  of  cases  then  is  m£.  For  example,  the  use  of  30  modes  and  3 
directions  gives  a  total  of  90  cases.  Rather  than  having  the  user  generate  the 
SULCASE  cards,  CASEGEN  generates  a  new  Case  Control  Data  Table  which  includes  the 
required  cards. 

The  DMAP  statement  for  CASEGEN  is 


CASEGEN  CAS ECC/ CASEDD /  C ,  ¥ , LMODES / V ,N,NDIR/V,N, NMODES  $ 


Parameter  LMODES  has  the  same  meaning  here  as  in  functional  module  GENPART  and 
must  appear  on  .  PARAM  bulk  data  card,  or  else  a  NASTRAN  fatal  error  will  result. 


NRLSUM 

Functional  module  NRLSUM  amputes  the  NRL  sum  stresses  and  forces  over  the  ra 
maximum  responses  for  a  given  direction  for  each  requested  element.  The  MRL  sum 
stress  for  a  given  stress  component  is 


j  i  max 


(ID 


where  S.  is  the  stress  component  for  the  1th  mode  and  S  is  the  maximum  of  these 
j  '  max 

stress  components  taken  ovet  all  modes  under  consideration.  The  Case  Control 
requests  for  stresses  and  forces  are  made  In  the  usual  wav,  except  that  S0RT2  format 
must  be  specified.  The  output  device  for  the  NRL  sums  (printer,  punch,  or  both)  will 
be  the  same  as  that  lor  the  standard  stresses  and  forces.  If  principal  stiesses  are 
computed  Cor  an  element,  they  will  be  computed  on  the  basis  of  the  NRL  sum  of  the 
normal  stresses.  For  the  EAR  element,  the  element  axial  stress  in  a  mode  will  be 
added  to  each  of  the  extenslonal  stresses  due  to  bending  in  that  mode.  The  NRL  sums 
will  then  be  computed  for  these  new  extensions',  stresses.  The  NRL  sums  corresponding 
to  the  printed  columns  headed  by  AXIAL  STRESS,  SA-MAX,  SB-MAX,  SA-HIN,  and  SB-NIN 
will  oe  set  to  0. 


In  seismic  analyses,  the  square  root  of  the  sums  of  the  squares  (SQRSS)  is  used 
rather  than  the  NRL  sums  for  the  stresses  and  forces.  The  user  may  select  the  latter 
method. 

The  DMAP  statement  for  NRLSUM  is 

NELSUM  0ES2 , 0EF2/NRLSTR, NRLFOR/V »N, NMODES/V , N, NDIR/C, Y , DIRECT  =  123/C,Y, SQRSS  =  0  $ 

Parameter  DIRECT  has  the  same  meaning  here  as  in  functional  module  GENCOS.  Integer 
parameter  SQRSS  indicates  whether  the  summing  process  should  use  NRL  sums  or  the 
SQRSS  method.  A  value  of  0  (the  default)  iudieates  NRL  sums;  a  value  of  1,  SQRSS. 

CQMBUGV 

COMBUGV  combines  the  direction-dependent  maximum  responses  in  a  number  of  ways. 

\ 

The  work  being  reported  is  -intended  for  DDAM  analyses,  but  seismic  analyses,  which 
make  use  of  a  similar  theory,  may  also  be  run.  In  seismic  analysis,  unlike  DDAM,  the 
maximum  responses  in  the  three  directions  for  each  mode  are  combined  into  one  total 
response  for  the  mode.  This  combination  may  be  performed  by  simply  adding  the 
absolute  values  of  the  maximum  component  responses  for  the  mode,  or  by  computing  the 
square  root  of  the  sum  of  the  squares  (SQRSS)  of  the  component  responses.  In  both 
cases,  the  result  is  a  matrix  in  which  each  column  represents  a  total  response  due 
to  a  mode.  These  responses  are  then  combined  by  the  SQRSS  method  over  the  modes  to 
give  a  final  response  vector.  Finally,  the  NRL  sums  of  the  displacements  are  also 
computed. 

The  DMAP  statement  for  COMBUGV  is 

COMBUGV  UGV/UGV ADD , UGV SQR , UGV ADC , UGVSQC , UGVNRL/ V , N , NMOD  ES / V , N , ND IR  $ 

Data  block  UGVADD  is  obtained  by  adding,  for  each  mode,  the  absolute  values  of  the 
component  responses  for  that  mode.  Data  block  UGVSQR  is  obtained  by  using  the  SQRSS 
method,  rather  than  by  adding.  Data  blocks  UGV ADC  and  UGVSQC  are  obtained  from 
UGVADD  and  UGVSQR,  respectively,  by  combining  the  total  modal  responses  using  the 
SQRSS  method.  Data  block  UGVNRL  contains  the  NRL  sums  of  the  displacements. 


II 


NASTRAN  INPUT  DETAILS 


A  complete  DDAM  analysis  with  NASTRAN  is  performed  in  one  normal  modes  analysis 
run  with  a  set  of  DMAP  ALTERS.  This  section  describes  the  input  details  for  such  a 
run. 

EXECUTIVE  CONTROL  DECK 

In  addition  to  standard  Executive  Control  Deck  cards,  the  Executive  Control 
Deck  for  the  normal  modes  analysis  must  include  the  proper  rigid  format  selection, 
SOL  3,0  and  the  following  DMAP  ALTER  package  (for  Level  17.5).  (The  numbers  to  the 
left  of  each  card  are  for  explanatory  purposes  only  and  are  not  actually  entered  on 
the  card.) 

1.  ALTER  133 

2.  GENCOS  BGPDT,  CSTM/DIRCOS/C,Y, SHOCK  =>  0/C,Y, DIRECT  -  123/ 
V,N,LUSET/V,N,NSCALE  $ 

3.  SAVE  NSCALE  $ 

4.  CHKPNT  DIRCOS  $ 

5.  DIAGONAL  MI/MID/C.N, SQUARE/ C,N,-1.  $ 

6.  MPYAD  MGG,  PHIG,/MP/C,N,0  $ 

7.  MPYAD  MP,  DIRCOS, /PMD/C,N,1  $ 

8.  MPYAD  MID,PMD,/PF/C,N,0  $ 

9.  CHKPNT  MID,  MP,  PMD,  PF  $ 

10.  DDAMAT  PF,  PMD/EFFW/C,Y,GG  «  386.4  $ 

11.  CHKPNT  EFFW  $ 

12.  LAMX,  , LAMA/LAMB/C, N,-l  $ 

13.  CHKPNT  LAMB  $ 

14.  GENPART  PF/RPLAMB,CPLAMB,RPPF,CPMP/C,Y ,LMODES/V,N,NMODES  $ 

15.  SAVE  NMODES  $ 

16.  PARTN  LAMB, CPLAMB , RPLAMB/ , , , OMEGA/ C,N,1  $ 

17.  CHKPNT  OMEGA  $ 

18.  PARAM  //C,N*GE/V,N, TEST/C, Y,LM0DES/V,N, NMODES  $ 

19.  COND  DDAM,  TEST  $ 

20.  PARTN  PF, ,RPPF/ ,PFR, ,/C,N,l  $ 


21.  CHKPNT  PFR  $ 

22.  EQUIV  PFR,PF  $ 

23.  CHKPNT  PF  $ 

24.  PARTN  EFFW,, RPPF/, EFFWR,, /C,N,1  $ 

25.  CHKPNT  EFFWR  $ 

26.  EQUIV  EFFWR, EFFW  $ 

27.  CHKPNT  EFFW  $ 

28.  PARTN  MP,CPMP,/,,MPR,/C,N,1  $ 

29.  CHKPNT  MPR  $ 

30.  EQUIV  MPR.MP  $ 

31.  CHKPNT  MP  $ 

32.  PARTN  PHIG,CPMP,/,,PHIGR,/C,N,1  $ 

33.  CHKPNT  PHIGR  $ 

34.  EQUIV  PHIGR, PHIG  $ 

35.  CHKPNT  PHIG  $ 

36.  LABEL  DDAM  $ 

37.  DESVEL  EFFW,0MEGA/SSDV,ACC,VWG,MINAC,MIN0W2/C,Y,GG=386.4/C,Y,VEL1/ 
C,Y,VEL2/C,Y,VEL3/C,Y, VELA/C, Y,VELB/C,Y,VELC/C,Y ,ACC1/ 
C,Y,ACC2/C,Y,ACC3/C,Y,ACCA/C,Y,ACCB/C,Y,ACCC/C,Y,ACCD  $ 

38.  CHKPNT  S  SDV , ACC , VWG , MINAC , MI N0W2  $ 

39.  DDAMAT  PF,MINAC/PVW/C,N,1.  $ 

40.  CHKPNT  PVW  $ 

41.  DDAMAT  PF,MIN0W2/PV0W/C,N,1.  $ 

42.  CHKPNT  PVOW  $ 

43.  DDAMPG  PHIG,PVOW/UGV/V,N,NMODES/V,N,NDIR  $ 

44.  SAVE  NMODES.NDIR  $ 

45.  CHKPNT  UGV  $ 

46.  DDAMPG  MP,PVW/PG/V,N,NMODES/V,N,NDIR  $ 

47.  CHKPNT  PG  $ 

48.  CASEGEN  CASECC/CASEDD/C,Y,LMODES/V,N,NDIR/V,N,NMODES  $ 

49.  EQUIV  CASEDD,CASECC  $ 

50.  SDR2  CASECC, CSTM, MPT, DIT,EQ EXIN, SIL, , ,BGPDT, ,QG,UGV,EST, ,/, 

OQG3,OUGV3,OES3*OEF3,/C,N,STATICS/V,N,NOSORT2  -  -1/C,N,-1  $ 
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51.  SAVE  N0S0RT2  $ 

52.  SDR3  0UGV3 , , 0QG3 , 0EF3 , 0ES3 , / 0UGV4 , , 0QG4 , 0EF4 , 0ES4 ,  $ 

53.  NRLSUM  0ES4 , 0EF4/NRLS  TR, NRLFOR/ V , N , NMODBS /V ,  N , NDIR/ C ,  Y , DIRECT  =  123/ 
C,Y,SQRSS  =  0  $ 

54.  OFP  NRLSTR, NRLFOR, , , ,  / /V,N, CARDNO  $ 

55.  SAVE  CARDNO  $ 

56.  COMBUGV  UGV/UGVADD,UGVSQR,UGVADC,UGVSQC,UGVNRL/V,N,NMODES/V,N,NDIR  $ 

57.  CASEGEN  CASECC/CASEEE/C.N, l/V,N,NDIR/V,N,NMODES  $ 

58 .  SDR2  CASEEE , CSTM, MPT , DIT , EQEXIN , S IL , , , BGPDT , ,QG, UGVNRL ,EST, ,/, , OUGV5 , , ,  / 
C,N,STATICS/V,N,NOSORT2/C,N,-l  $ 

59.  SAVE  NOSORT2  $ 

60.  OFP  0UGV5,,,,,//V,N, CARDNO  $ 

61.  SAVE  CARDNO  $ 

62.  ENDALTER 

The  following  notes  on  the  functional  module  executions  are  keyed  to  the  cards  by 
corresponding  number. 

2.  Computes  direction  cosine  matrix  [V]  in  Equation  (la). 

5.  Creates  a  diagonal  matrix,  consisting  of  the  diagonal  of  the  modal  mass 
matrix,  and  inverts  it.  The  new  matrix  is  used  in  Equation  (2). 

6.  Computes  for  later  use. 

T 

7.  Computes  [ <t> ]  [M]  [V]  as  described  in  Equation  (la). 

8.  Computes  matrix  of  participation  factors  [P  ]  (Equation  (2)). 

Si 

10.  Computes  effective  masses  and  weights  in  Equation  (4). 

12.  Creates  a  matrix  of  the  information  in  the  Real  Eigenvalue  Table  for 
later  use  in  Equation  (5) . 

14.  Creates  partitioning  vectors  which  will  be  used  to  create  a  vector  of 
natural  circular  frequencies  from  a  matrix  of  miscellaneous  eigenvalue 
results.  Additionally,  if  the  number  of  modes  to  be  used  in  computing 
maximum  responses  is  less  than  the  number  computed  in  the  normal  modes 
analysis,  other  partitioning  vectors  are  created  to  reduce  the  orders 
of  a  number  of  matrices. 

16.  Creates  the  vector  of  natural  circular  frequencies. 

18.  Compares  the  number  of  desired  modes  (LMODES)  and  the  number  of  computed 
modes  (NMODES). 
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19.  If  LMODES  >  NMODES,  skips  to  36. 

20-35.  Reduces  orders  of  several  matrices  from  NMODES  to  LMODES. 

37.  Computes  shock  spectrum  design  velocities  and  accelerations,  as  given  in 
Equations  (6)  through  (8) .  In  addition,  matrices  corresponding  to 
Equations  (9)  and  (9a)  are  created  for  use  in  Equations  (5)  and  (5f), 
respectively.  If  Equations  (6)  through  (9)  do  not  represent  the  desired 
forms  for  velocities  and  accelerations,  matrix  MINAC  or  MIN0W2  may  be 
directly  specified  on  DMI  bulk  data  cards  and  functional  module  DESVEL  may 
be  deleted.  MINAC  and  MIN0W2  must  be  of  order  LMODES  x  L;  LMODES  is 
explained  in  6  above,  and  L  is  the  number  of  desired  directions. 

39.  Creates  the  outer  product  of  Equation  (10). 

41.  Creates  the  matrix  of  Equation  (10a). 

43.  Computes  the  LMODES  *  L  matrix  of  direction-dependent  maximum  responses. 

46.  Creates  the  LMODES  *  L  static  load  matrix  as  in  Equation  (5). 

48.  Generates  a  new  Case  Control  Data  Table  which  includes  the  (LMODES  *  L) 
subcases. 

50.  Computes  stresses  due  to  each  maximum  response. 

52.  Converts  stresses  in  5C.  from  S0RT1  to  SORT2. 

53.  Computes  the  NRL  sum  or  SQRSS  stresses  and  forces  for  each  requested 
element . 

54.  Outputs  the  NRL  sum  stresses  and  forces  to  the  printer  and/or  punch,  as 
requested  in  the  Case  Control  Deck. 

56.  Computes  various  combinations  of  the  component  maximum  responses. 

57-61.  Prepares  and  prints  file  of  NRL  sum  displacements. 

CASE  CONTROL  DECK 

Although  the  usual  selections  may  be  made,  two  requirements  are  imposed: 

1.  No  subcases  are  to  be  specified, 

2.  Stress  and  force  selections  must  request  S0RT2  format. 

This  last  requirement  will  force  all  output  selections,  e.g.,  displacements,  applied 
loads,  etc.,  to  be  in  SORT2  format.  Also,  the  NRL  sum  stresses  and  forces  will  be 
printed  and/or  punched,  as  requested  in  the  corresponding  STRESS  and  FORCE  requests. 
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BULK  DATA  DECK 

The  values  of  a  number  of  parameters  special  to  DDAM  must  be  specified.  For 
those  parameters  with  no  default  values  and  for  parameters  for  which  the  default 
values  are  to  be  overridden,  PARAM  bulk  data  cards  will  be  required.  The  parameters 
are  as  follows: 

1.  SHOCK  -  The  nonnegative  integer  value  of  this  parameter  is  the  identifica¬ 
tion  number  of  the  coordinate  system  which  defines  the  shock  directions.  A  non-zero 
value  requires  definition  of  the  system  on  a  CORDij  card.  A  zero  value  implies  the 
basic  coordinate  system  with  shock  directions  X,  Y,  and  Z.  The  default  value  is 
zero.  The  value  of  SHOCK  should,  in  most  cases,  be  the  same  as  the  displacement 
coordinate  system  identification  number  for  the  grid  points. 

2.  DIRECT  -  This  parameter  may  have  one  of  the  following  integer  values:  1,  2, 
3,  12,  13,  23,  or  123.  The  default  value  is  123.  The  value  of  DIRECT  indicates 
which  directions  of  coordinate  system  SHOCK  are  to  be  considered.  For  example,  if 
DIRECT  *»  123,  then  all  three  directions  will  be  used.  If  DIRECT  =  13,  only  two 
directions  will  be  used,  the  first  and  the  third. 

3.  GG  -  This  parameter  is  the  acceleration  due  to  gravity.  The  default  value 
is  386.4. 

4.  LMODES  -  The  integer  value  of  this  parameter  is  the  number  of  lowest  modes 
to  be  used  in  the  static  analyses.  This  number  may  be  less  than  the  number  of  modes 
computed  in  the  normal  modes  analysis.  No  default  value  is  provided  so  the  value  of 
this  parameter  must  be  given  on  a  PARAM  bulk  data  card  or  else  a  fatal  message  will 
result . 

5 .  VEL1 , VEL2 , VEL3 , VET  .A , VELB , VELC , ACC1 , ACC2 , ACC 3 , ACCA , ACCB , ACCC , ACCD  - 

The  real  values  of  these  parameters  control  the  computation  of  the  shock  spectrum 
design  values  for  velocity  and  acceleration.  These  parameters  are  defined  by 
Equations  (6)  through  (8)  and  further  explained  on  pages  7  and  8.  No  default  values 
for  any  of  these  parameters  are  provided,  so  a  PARAM  bulk  data  card  for  each 
parameter  must  be  included  in  the  Bulk  Data  Deck. 
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